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A formula to investigate wave effect in multi-lens system is presented on the basis of path integral 
formalism by generalizing the work by Nakamura and Deguchi (1999). Wave effect of a system with 
two lenses is investigated in an analytic way as a simple application to demonstrate usefulness of 
the formula and variety of wave effect in multi-lens system. 
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I. INTRODUCTION 
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Recently several authors have investigated wave effect in gravitational lensing phenomenon, motivated by the 
possibility that the wave effect might be detected in future gravitational wave observatories 

wave effect can be investigated by solving wave equation, in principle. However, it is difficult to solve the wave 
Oh- equation analytically for general configuration of lenses, excepting the simple case of a single Schwarzschild lens (e.g., 
see also and references therein, cf. Q). Indeed, investigation of the wave effect so far is restricted to special 
case of single lens model with the spherically symmetry, i.e., the Schwarzschild lens model and the singular isothermal 
^""J | sphere lens model (see Q). 

Nakamura and Deguchi developed an elegant formalism for gravitational lens using the path integral approach 
. The primary purpose of the present paper is to derive a useful formula to investigate the wave effect in general 
. multi-lens system, by generalizing the formalism by Nakamura and Deguchi. We also apply it to a system with two 
\^ ' Schwarzschild lenses to demonstrate usefulness of the formula. This paper is organized as follows: In §. 2, we present 
Q\ , a generalized formula for multi-lens system. In §. 3, an application of the formula to a simple configuration with two 
• lenses is considered. §. 4 is devoted to summary and conclusion. We use the convention c = 1. 
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We start by reviewing the path integral formalism for gravitational lens phenomenon Q . We consider the Newtonian 
II i spacetime with the metric 

^ ■ ds 2 = -(1 + 2U(r, 9, Lp))dt 2 + (1 - 2U(r, 9, ip)){dr 2 + r 2 d6 2 + r 2 sin 2 9dip 2 ), (1) 

where U(r, 9, ip) is the Newtonian potential. Propagation of massless field <j> is described by the wave equation 
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II. GENERALIZED FORMULATION 



where g is the determinant of the metric and r represents the spatial coordinates (r, 9, tp). We consider a monochro- 
- ' matic wave from a point source with the wave number k. We set 

^,t) = ^(r)e-'^, (3) 

where A is a constant, then Eq. yields 

Assuming that the first term is negligible compared with the second term and 9 <C 1, Eq. reduces to 

■2kU(r,9,tp)F. (5) 

Due to analogy of Eq. JSJ with the Schrodinger equation, using the path integral formulation of quantum mechanics, 
the solution can be written as follows, 
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F{r,9,ip)= JvOexp ik dr Qr 2 9 2 (r) - 2U(r,Q{r)) 
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where the dot denotes the differentiation with respect to r and O is used to represent the variables (9, ip), which are 
related by 9 = (8 cos ip, 8 sin ip). The expression © means the sum of all possible path 0(r), fixing the initial point 
(source's position) and the final point (observer's position). 

Let us consider multi-lens system, including n lenses as shown in Fig. 1, in which the source is located at the origin 
of coordinate and an observer is located at (r, 0) — (r/v,0jv). The radius between the source and the observer is 
discretized by N segments, and each discrete radius from the source is labeled by rj with j from 1 to N. Separation 
neighboring two discrete radii is e. We assume that the n lenses are located at the radius r = n m with m from 1 to n, 
and that the thin lens approximation is valid for the lenses. In this case the explicit expression of Eq. (JSJ) is written as 
Eq. I|A2|) in Appendix (see also [6J). After some computation, the path integral expression reduces to (see Appendix 
for details) 
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with 



^(0,J = 2 f ^ + drU(r,B), 



(8) 



where Ti m _ lt i m = ri m — rj m _ l , n„ +1 = rjv, and 0; m is the variable on the mth lens plane. As expected, stationary 
condition of the phase of the integrant in Eq. Q reproduces the lens equation in the multidens system 
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for each m! = 1, • • ■ , n, which hints at a way to the geometrical optics limit. The Gaussian approximation around a 
stationary solution yields the result in the geometrical optics limit (see also |(|). 

III. APPLICATION TO A SIMPLE CONFIGURATION 

In this section we consider a simple case with two lenses. Fig. 2 shows the configuration: The Schwarzschild lenses 
with mass Mi and Mi are located at the radius r\ and r%, respectively. The model considered here is not general, 
because the source and the two lenses are arranged to be on a straight line. However, this simplification allows us to 
perform integration of Eq. Q analytically. We start by rewriting Eq. @ 
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with 



^(e j ) = 4GM i ln(|6 J -|) : 



for j = 1, 2, where rij = rj — rj and the possition of an observer is (r3, 03). From Eq. Hl()|l we have 
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where \Qj\ = 9j and Jo{y) is the Bessel function of the first kind. Integration with respect to 9\ can be performed 
(see E3) 
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where a is a real number which represents a constant phase and i-Fi(a, 6 ; y) is the Kummer's function. With the use 
of the definition of the Bessel function 



we have (see [TT| 
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where we defined 



x {xz) L 2 F 1 (1 - 2ikGM 1 , 1 + L- 2ikGM 2 , 1 ; 1 - z) , (15) 
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a' is a real constant and 2 Fi(a, b,c ;y) is the Hypergeometric function. In the limit 9 3 = (x — 0), Eq. (fT3f) reduces 
to 

F = e la ' e' nkG{Ml+M2) T{l-2ikGM 1 )T{l~2ikGM 2 ) 

xza^i (1 - 2ikGM u 1 - 2ikGM 2 , 1 ; 1 - z) , (18) 

and we have 

If - x x ^CM, * \*K (1 - ^G Mll 1 - 2, fc GM 2 , 1 ; 1 - ztf . (19) 

We consider the coincidence limit that the distance between the two lenses becomes zero, i.e., r x = r 2 . In this case, 
from previous investigation (e.g., [T(ijp. i* 1 should be 

T(x) = e - fc G(Af 1+ Af 2 ) r(1 _ 2ikGMx)xFi (1 - 2ifcG(Mi + M 2 ), 1 ; -is) , (20) 

excepting a constant phase factor. We compare our result with Eq. (|20(1 . First, let us consider the limit r\ = r 2 , i.e., 
z = of Eq. (|15f) . Using the mathematical formula 

T(l - 2ikGM 1 )T(l ~ 2ikGM 2 )z 2 F 1 (1 - 2ikGM x , 1 - 2ikGM 2 , 1 ; 1 - z) 

= z 2ifeG(M 1+ M 2 ) r ^ _ 2ikG(M 1 + M 2 )) 

x 2 ^i (2ikGM l7 2ikGM 2 , -1 + 2ikG(M 1 + M 2 ) ; z) 

T(l - 2i/cGMi)r(l - 2ifcGM 2 )r(-l + 2ikG(M 1 + M 2 )) 
+Z T(2ikGM 1 )T(2ikGM 2 ) 

x 2 ^i (1 - 2ikGM u 1 - 2ikGM 2 , 2 - 2ikG(M 1 + M 2 ) ; z) , (21) 
in the case 83 = (x — 0), we can easily have 
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This is the expected result. Note that |F| 2 is regarded as the magnification factor. Fig. 3 plots R — |F/JF(0)|, as a 
function of z, where we fixed kGMi = kGMi = 1. Thus the maximum magnification depends on lens configuration 
significantly. Fig. 4 plots R = \F/F(x)\ as a function of x when fixing z = 0.5 and kGM\ = kGMi = 1, which indicates 
that interference pattern also depends on lens configuration [l^. In these figures we have fixed kGM\ — kGMi = 1, 
but result depends significantly on the parameters kGMi, kGMi and z, as shown in Fig. 4, which suggests variety of 
wave effect depending on lens configuration. 



IV. SUMMARY AND CONCLUSION 



In the present paper, we have presented a general formula to investigate wave effect in multi-lens system, which has 
been derived with the use of path integral approach. The formula is expressed in terms of integration with respect to 
variables of lens planes. It is difficult to perform the integration in general cases, but a system with two Schwarzschild 
lenses is an example for which the integration can be performed in an analytic way. The model of considered in the 
present paper is simplified and limited, however, it suggests variety of wave effect in multi-lens phenomenon. It is 
required to develop a numerical method to perform integration in general lens configuration in future work. 
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APPENDIX A: DERIVATION OF EQ. © 



In this Appendix, we review derivation of Eq. J7J from Eq. JHJ). We consider the configuration depicted as Fig. 1. 
The source is located at the origin of coordinate and the position of an observer is specified by (r, 9) = (rjv,9jv)- 
The space between the source and the observer is discritized by N — 1 planes. The radius of each plane is labeled by 
rj for j = 1, • • • , N — 1. Note that specifies the plane of the observer. 9j is (angle) variable on the jth plane. We 



consider the system with n lenses, and assume that mth lens is located at the radius 
the thin lens approximation, we introduce the two dimensional potential 



Assuming the validity of 



^(e,J = 2 / ' + drU(r,&), 
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for m = 1, • • * , n, respectively. In this case, the path integral formula JSJ) can be written as 
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where the normalization is chosen 
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2nie 



kr 



(A3) 



so that F — 1 in the limit "0 = 0. Then, Eq. l|A2j) is rephrased as 
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where e is the separation between two neighboring planes. With the use of the following equality, which can be proven 
by the mathematical induction, 
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which is equivalent to Eq. Q 
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FIG. 1: Configuration of multi-lens system and coordinates for path integral formula. 
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FIG. 2: Lensing system considered in section 3. 
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FIG. 3: R = \F/T(0)\ as a function of z with fixing fcGMi = kGM 2 = 1. 




FIG. 4: R = \F/T{x)\ as a function of x with fixing z = 0.5 and fcGMi = kGM 2 = 1. 



